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DECOUPLING  ION  EFFICIENCY  FROM  IMPEDANCE  IN 
PULSED  POWER  DIODES 


I.  INTRODUCTION 

Nearly  a  decade  has  elapsed  since  the  Intense  Ion  Beam  Research  Group  at 
NRL  first  clearly  outlined  the  Important  role  light  ions  play  in  the  physics 

of  pulsed-power  diode  operation.^  Building  upon  a  theoretical  understanding 

2 

of  the  nature  of  electron  flow  in  the  diodes  coupled  with  a  storehouse  of 
experimental  observations'5 »  the  group  developed  a  semi -empirical  prediction 
of  net  diode  current  in  the  form  of^ 


I  +  I. 
e  i 
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where  y  ■  1  +  ,  R  ■  cathode  radius,  D  *  axial  anode-cathode  (A-K)  gap, 
V  »  diode  voltage,  and  m^  »  ion  mass.  Implicit  in  this  formula  is  an  ion-to- 
electron  current  ratio  which  may  be  approximated  by6 


I  c  D 

e 


(2) 


where  v^  is  the  maximum  ion  velocity  and  the  electrons  are  assumed  to  have 
been  Instantaneously  accelerated  to  near-light  velocities.  According  to  these 
two  expressions  and  for  a  fixed  machine  voltage,  the  diode  impedance,  V/I, 
will  vary  inversely  with  the  aspect  ratio,  R/D,  while  the  relative  proportions 
of  diode  current  carried  by  the  ions  will  vary  directly  with  R/D. 

In  order  to  understand  the  physics  of  Eqs.  (1)  and  (2)  it  is  useful  to 

refer  to  a  schematic  picture  of  a  representative  axial  geometry  pulsed-power 

diode.  (The  distinction  between  axial  and  radial  geometries  is  given  by  the 

orientation  of  their  respective  primary  anode-cathode  gap.  The  radial 

geometry  will  be  discussed  later.)  Figure  1  depicts  the  nature  of  electron 
Manuscript  approved  July  22, 1983.  1 


and  ion  flow  in  the  diode.  First  consider  the  case  of  electrons  flowing  in 
the  absence  of  ions,  1^  -  0.  It  has  been  shown ^  that  the  largest  fraction  of 
electron  current  is  emitted  from  the  "edge"  portions  (i.e.,  r  -  R)  of  the 
cathode.  The  emitting  surfaces  at  those  large  radii  benefit  from  two- 
dimensional  enhancement  of  the  normal  electric  field  while  the  regions  at 
lower  radii  experience  field  strength  reductions  due  to  the  large  amounts  of 
electron  space  charge  emitted  at  larger  radii.  Thus,  the  combination  of 
enhanced  electron  emission  current  density  with  the  larger  emission  area  makes 
the  behavior  of  the  electron  flow  at  r  ■  R  crucial  to  the  final  operating 
state  of  the  diode  as  a  whole.  Ihe  azimuthal  magnetic  field  generated  by  the 
diode  current  flowing  within  that  radius  plays  a  role  equal  to  that  of  the 
electric  field  in  determining  the  trajectories  of  those  outer  electrons.  It 
is  this  feedback  loop  between  the  diode  current,  electron  emission 
at  r  *  R,  and  BQ(r*R)  that  provides  the  rational  for  Eq.  (1).  To  be  specific, 
if  the  net  current  exceeds  some  "critical"  value,  Icrit ,  then  BQ  in  the  outer 
radial  regions  will  cause  the  mean  Larmor  radius  of  electrons  there  to  be  less 
than  the  anode-cathode  gap  spacing,  D.  In  that  case,  those  outer  electrons 
will  not  be  able  to  cross  directly  from  the  cathode  to  the  anode  but  rather 
will  be  turned  around  toward  the  cathode  again  at  some  lower  radius  where 
their  additional  negative  space  charge  will  cause  reduced  electron  emission. 
The  reduced  electron  current,  Ig,  will  lower  Bg(r«"R)  thus  enlarging  the  mean 
Larmor  radius  until  it  again  exceeds  the  gap,  D,  and  allows  direct  electron 
streaming  to  the  anode.  Further  weakening  of  Bg  beyond  that  value  is  resisted 
by  the  again  unimpeded  electron  emission  at  lower  radii.  This  reasoning 
argues  in  favor  of  an  equilibrium  I@  value  in  the  neighborhood  of  an  Icr^t  for 
which  the  resultant  mean  electron  Imrmor  radius  at  r»R  just  equals  the  gap 

Q 

size,  D.  After  some  algebra  this  "critical  current"  in  kA  may  be  written  as 


where  and  are  the  standard  relativistic  quantities  for  electrons  in  the 

Q 

diode  at  voltage,  V.  Although  recent  theoretical  and  experimental  work 
indicates  that  attention  to  electrode  shaping  and  to  ion  space  charge 
distribution  can  boost  electron  currents  almost  a  factor  of  two  above  lcrit» 
it  still  serves  as  a  scaling  guideline  for  diode  operating  parameters.  Thus 
the  premise  for  Eq.  (1)  remains  intact  aside  from  a  possible  change  in  the 
proportionality  factor. 

The  rationale  underlying  Eq.  (2)  can  likewise  be  explained  by  reference 
to  Fig.  1.  Stated  simply,  the  relative  emission  of  ions  and  electrons  from 
their  respective  electrodes  is  a  function  of  the  electric  field  strengths 
normal  to  the  electrode  surfaces  in  the  A-K  (anode-cathode)  gap.  Those 
fields,  in  turn,  depend  not  only  upon  the  diode  voltage,  V,  and  A-K  gap 
spacing,  D,  but  also  upon  the  relative  amounts  of  ion  and  electron  space 
charge  in  the  gap.  Therefore,  the  ion-to-electron  current  ratio  will  depend 
inversely  on  the  relative  "dwell  times"  of  the  respective  species  in  the 
active  A-K  gap  region.  The  dwell  time  ratio  may  be  expressed  as  Tg/T^.  In 
the  pinched-beam  mode,  the  predominant  flow  of  electrons  from  the  outer 
cathode  radii  will  travel  about  a  distance,  R,  through  the  A-K  gap  at  about  a 
uniform  velocity,  c,  for  diode  voltages  of  a  megavolt  or  more.  The  emitted 
ions,  on  the  other  hand,  are  hardly  perturbed  at  all  from  simple  paraxial, 
straight-line  flow  through  the  gap  across  its  width,  D.  For  a  diode  potential 
difference,  V,  they  will  attain  velocities  of  [2  eV/m^]  ^ .  This  reasoning 
leads  to  a  restatement  of  Eq.  (2)  in  the  form 


The  situation  changes  significantly  when,  instead  of  considering  the 

axial  diode  configuration  of  Fig.  1,  one  turns  to  a  radial  diode  configuration 

such  as  that  in  Fig.  2.  In  this  new  orientation,  the  total  diode  current  can 

be  assumed  to  ideally  split  evenly  on  either  side  of  the  diode  centerplane. 

(The  stability  of  such  an  even  current  division  remains  to  be  rigorously 

tested.)  These  currents  will  generate  an  azimuthal  magnetic  field,  B^,  in  the 

A-K  gap  analogous  to  that  present  in  the  axial  diode.  Urns,  it  can  be 

expected  that  the  current  on  either  side  of  the  half-plane  will  face  the  same 

Icr£t  limitation  as  that  expressed  in  Eq.  (3)  for  the  total  axial  diode 

current.  Therefore,  in  the  case  of  the  radial  diode,  I..  ,  "I  . .  and  the 

*  '  diode  crit 

diode  impedance  is  still  inversely  linked  to  the  ratio  of  R/D. 

This  correspondence  between  the  two  diode  geometries  breaks  down  when  one 

now  attempts  to  approximate  the  ion-to-electron  current  ratio.  It  can  again 

be  assumed  that  I^/I^  will  be  inversely  proportional  to  the  respective  species 

"dwell  times".  Here  again  the  ion  dwell  time,  t^,  can  be  approximated 

Vo 

by  [eV/(2m^)]  z/ D.  The  electron  flow,  on  the  other  hand,  is  no  longer 
collapsing  radially  inward  under  the  action  of  and  Ez.  Rather,  the 
electric  field  is  now  dominated  by  Er  so  that  the  electrons  are  now  being 
forced  to  drift  to  the  centerplane.  They  must  therefore  cover  a  distance,  L, 
Instead  of  R  and  they  now  enjoy  dwell  times  more  closely  approximated 
by  ■  L/c.  Thus,  the  ion  generation  efficiency  of  the  diode  no  longer 
depends  upon  R/D  and,  thereby,  is  no  longer  directly  linked  to  the  diode 
Impedance.  The  implications  of  such  a  decoupling  are  highly  significant  to 
researchers  seeking  to  maximize  either  ion  beam  or  electron  beam  production 


for  a  diode  with  fixed  impedance  parameters. 

This  report  will  present  numerical  simulations  of  example  diode 
geometries  for  both  the  enhanced  ion  current  as  well  as  for  the  enhanced 
electron  current  cases.  The  numerical  findings  will  be  compared  to 
theoretical  predictions  as  well  as  to  the  experimental  observations  currently 
available.  The  geometry  for  the  enhanced  ion  current  case  was  based  on  a 
design  proposed  for  the  joint  NRL-Sandia  SHEPERD  (Self -magnetically  Healing 
Equatorial  Pinch  Electron-ref lexing  Radial  Diode)  project.^  The  simulations 
confirmed  assumptions  regarding  the  nature  of  ion  and  electron  flow  in  the  A-K 
gap.  Furthermore,  experiments  conducted  on  the  Gamble  II  pulsed  power 
generator  to  test  the  geometry  confirmed  the  strong  dependence  of  the 
impedance,  Z,  on  L.  The  second  set  of  simulations  focused  on  the  problem  of 
minimizing  ion  production  and  maximizing  the  electron  current  in  a  low 
impedance  diode.  Here,  the  simulations  demonstrated  unequivocally 
that  I^/Ie  was  tracking  L/D  and  not  R/D,  For  this  enhanced  electron  case,  a 
subset  of  simulations  was  conducted  to  demonstrate  that  the  physics  of  the  L/D 
dependence  could  be  translated  to  a  diode  in  axial  geometry  as  well  by 
properly  routing  the  diode  return  current  flow  path  with  respect  to  the  A-K 
gap.  Implications  of  the  overall  results  for  several  applications  are  then 
discussed. 

II.  ENHANCED  ION  CURRENT 

Much  research  in  the  field  of  light  ion  inertial  confinement  fusion  has 

been  devoted  to  the  problem  of  boosting  ion  efficiency  in  high  impedance 
11  12 

diodes.  A  switch  from  axial  to  radial  geometry  holds  promise  for  being 

the  final  solution. 

To  test  the  physical  assumptions  regarding  the  nature  of  electron  flow  in 


Che  radial,  pinched-beam  model,  the  geometry  depicted  in  Fig.  3  was  chosen  as 

a  subject  for  numerical  simulation.  The  cathode  radius,  R,  was  taken  as  5.0 

cm  to  equal  the  gap  half-width,  L,  of  5.0  cm.  The  radial  A-K  gap  itseir  was 

fixed  at  0.2  cm  while  the  side-gap  in  the  magnetically  insulated  triplate 

power  feed  was  set  at  0.9  cm.  Two  types  of  anode  surfaces  were  studied.  The 

first  was  a  completely  solid,  electron-absorbing  surface  spaced  0.2  cm  away 

from  the  cathode  and  along  which  the  axial  extent  of  the  ion  emission  could  be 

varied.  (Note  that  complete  mirror  symmetry  about  the  diode  center  plane  was 

assumed  at  all  times.)  The  simulation  of  the  diode  was  then  repeated  assuming 

this  time  that  the  front  anode  surface  was  a  thin  (~  5  mil  polyethylene)  foil 

held  at  full  anode  potential  and  through  which  1.5  MeV  electrons  could  easily 

travel.  This  foil  becomes  a  plasma  3heet  when  the  diode  is  in  operation. 

Behind  the  foil  was  a  full  0.4  cm  radial  gap  separating  it  from  the  solid 

anode  surface.  This  is  the  configuration  shown  in  Fig.  3.  There  was  assumed 

to  be  zero  electric  field  in  the  gap  behind  the  foil  but  the  full  B  due  to 

any  axial  currents  below  that  radius  could  be  felt  there.  Thus,  electrons 

penetrating  through  the  foil  from  the  A-K  gap  could  be  expected  to  experience 

a  bending  force  due  to  the  BQ  behind  the  foil.  If  that  magnetic  field  is 

strong  enough,  the  electrons  could  be  completely  turned  around  and  reinjected 

through  the  foil  back  into  the  radial  A-K  gap.  This  is  completely  analogous 

to  the  electron  behavior  in  the  axial  pinch-reflex  diodes  pioneered  by 
13 

NRL.  Such  electron  reflexing,  if  shown  to  occur  for  this  geometry,  can  be 
expected  to  further  enhance  the  production  of  light  ion  flux  from  the  diode. 

To  examine  the  behavior  of  this  new  diode  configuration  for  the  cases 
described  above,  numerical  simulations  were  conducted  using  NRL's  PREMAS 
computer  code.  As  the  successor  to  DI0DE2D,  the  new  PREMAS  code  embodies  all 
of  its  .'edeces'  s  accuracy  and  efficiency  in  characterizing  steady  state. 


self-consistent  ion  and  electron  flows  in  pulsed  power  diodes.  In  addition  to 
those  attributes,  it  also  has  an  unparalleled  flexibility  in  dealing  with  a 
huge  array  of  axial  and  radial  diode  geometries.  Finally,  it  boasts  a 
completely  new  diagnostics  package  based  on  the  powerful  DISSPLA  graphics 
software  system.  The  details  of  the  code's  field-solving  and  particle-pushing 
algorithms  may  be  found  elsewhere.*^ 

It  is  sufficient  here  to  point  out  that  PREMAS  calculates  equilibrium 
electric  and  magnetic  field  strengths  over  an  NZ  x  MR  mesh  of  discrete  data 
points  on  a  predetermined  computational  region  lying  on  a  arbitrary  R  -  Z 
planar  cross  section  passing  through  the  diode's  centerline.  (Mote,  however, 
that  the  lower  boundary  of  the  computational  region  itself  need  not  correspond 
to  the  centerline.  Rather,  it  is  free  to  lie  at  some  nonzero  value  of  the 
physical  diode  radius.)  Complete  azimuthal  symmetry  is  assumed.  A  finite 
number  of  macroelectrons  and  macroprotons  having  correct,  physical  charge-to- 
mass  ratios  are  advanced  times tep-by -times tep  across  the  mesh  in  a 
relativistically  covariant  manner.  A  self-consistent,  steady-state  solution 
is  sought  both  for  field  structures  as  well  as  for  particle  flows.  The  time- 
dependent  Maxwell's  Equations  are  not  observed. 

For  this  specific  problem,  the  physical  volume  simulated  by  the  PREMAS 
code  corresponded  to  the  "left  half"  of  the  A-K  gap  region  of  the  Fig.  3 
diode.  This  numerical  model  region  is  depicted  in  Fig.  4.  Numerical  cell 
spacings  of  ARxAZ  *  0.01  x  0.023  cm  were  chosen  to  fill  the  computational 
region.  Thus  the  number  of  axial  and  radial  grid  points  are  256  and  60, 
respectively,  to  fill  the  5.9  by  0.6  cm  physical  space.  Allowing  a  monolayer 
of  "quard  cells"  surrounding  the  entire  perimeter  of  the  region  gives  a  total 
cell  count  of  NZ2  x  NR2  ■  258  x  62  ■  15,996  cells.  In  all  cases,  the  lower 
boundary  was  the  cylindrical,  solid  cathode  surface  at  R  *  5.0  cm  which  was 


held  fixed  at  zero  voltage.  Similarly,  the  "left"  boundary  was  always  the 
planar,  solid  cathode  surface  designated  to  be  z  =  0.0  and  also  held  to  zero 
voltage.  For  the  solid  anode  case,  the  upper  boundary  was  at  R  *  5.2  cm  while 
for  the  case  with  the  reflexing  foil  it  was  set  back,  to  R  *  5.6  cm.  For  both 
cases  the  electrostatic  boundary  potential  there  was  graded  linearly  for  Z  = 
0.0  to  Z  *  0.9  cm  from  $  =  0  to  $  =  +  1.5  MV.  The  rest  of  the  upper  boundary 
was  maintained  at  full  anode  potential  at  *  +1.5  W.  Ihe  treatment  of  the 
"right"  boundary  lying  on  the  diode  centerplane  is  a  bit  more  complicated.  As 
already  stated,  the  right  half  of  the  diode  is  assumed  to  be  a  mirror  image  of 
the  left  half  across  the  centerplane.  Hius,  from  a  computational  standpoint 
this  "right"  boundary  is  perfectly  reflecting  to  particle  flows.  The  presence 
of  space  charge  and  current  there  also  prohibits  potential  values  from  being 
fixed  along  it.  The  value  of  <f>  will  change  there  during  the  simulation  so  it 
is  allowed  to  "float"  from  timestep  to  timestep  in  a  slightly  damped 
fashion.  This  "floating  boundary  feature  is  enforced  only  in  the  A-K  gap 
itself  (i.e.,  between  R=«5.0  and  R*5.2)  for  both  cases.  For  the  anode  foil 
case,  of  course,  the  "right"  boundary  between  R=5.2  and  R=5.6  runs  inside  the 
anode  conductor  and  is  therefore  maintained  at  a  constant  <j>  *  +1.5  W. 

The  first  case  involved  an  analysis  of  diode  operation  with  a  solid, 
electron-absorbing  anode  surface  in  place  of  the  foil.  Ion  emission  along 
that  surface  was  turned  on  in  successive  stages  to  examine  the  steady  state 
electron  flow  for  each  stage.  Initially,  only  electron  emission  was 
permitted.  The  equilibrium  flow  is  depicted  in  Fig.  5a.  (Note,  again,  that 
only  half  of  the  diode  is  treated  numerically  and  that  mirror  symmetry  about 
the  centerplane  is  assumed.)  In  this  figure,  as  well  as  in  the  others  that 
follow,  the  radial  dimension  scale  has  been  exaggerated  to  show  maximum  detail 
in  the  flow  pattern.  Two  significant  results  are  observed.  First,  the  net 
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electron  current  in  this  half-diode  is  0.85  MA  which  is  almost  identical  to 
the  critical  current  of  0.81  MA.  In  addition,  there  is  a  slight  axial 
deflection  of  the  electron  sheath  as  it  crosses  the  A-K  gap.  Ihe  net  Inward 
deflection  measured  0.21  cm  -  almost  identical  to  the  gap  width  of  0.20  cm. 
This  amounts  to  a  near-perfect  quarter-cycle  of  electron  gyromotion. 

The  next  question  to  be  considered  concerned  the  possible  effects  of 
limited  ion  emission  along  the  anode  surface.^  Suppose  the  electron  beam, 
through  one  physical  mechanism  or  another,  creates  a  plasma  along  the  anode 
surface  it  is  striking.  To  simulate  that  effect,  proton  emission  was  "turned 
on"  along  a  0.34  cm  wide  strip  from  Z*1.26  to  Z»1.60  cm  on  the  solid  anode 
surface  immediately  surrounding  the  electron  impact  zone.  The  net  result  is 
shown  in  Fig.  5b.  Thanks  to  even  the  low  ion  current  of  75  kA,  enough 
positive  space  charge  has  been  introduced  into  the  A-K  gap  near  the  electron 
emission  strip  to  boost  Ig  to  1.075  MA.  This  brings  the  total  half -diode 
current  to  over  1.4  x  lcrit*  As  would  be  expected  from  this  increased  diode 
current,  the  axial  pinching  of  the  electron  sheath  is  more  pronounced.  The 
net  displacement  is  now  1.42  cm  toward  the  centerplane.  In  the  next  step,  the 
width  of  the  ion  emission  strip  was  nearly  tripled  so  that  it  extended  from 
Z*1.26  to  Z-2.30  cm.  The  change,  as  indicated  in  Fig.  5c  was  minor.  The  net 
ion  current  did  not  quite  double,  but  rather  rose  to  140  kA.  Furthermore,  in 
response  to  the  increased  strength  of  B  along  the  outer  cathode  surface,  the 
net  electron  current  actually  dropped  to  1.03  MA.  Still,  the  total  half -diode 
current  remained  at  over  1.4  x  Icr^t .  It  was  clear  from  this  new  steady-state 
finding  that  the  electron  beam  would  not  simply  march  axially  inwards  in-step 
with  an  expanding  zone  of  ion  emission.  Some  minimum  threshold  for  beam 
pinching  must  first  be  surpassed.  This  threshold  is  clearly  more  than  1.4  x 


In  order  Co  pass  Che  pinching  threshold,  ion  emission  was  curned  on  along 
Che  entire  solid  anode  surface  from  Z-1.26  inward  to  the  diode  centerplane. 

The  change  was  dramatic.  As  shown  in  Fig.  5d,  Che  electron  sheath  displayed 
classic  ExB  cycloids  while  converging  to  the  diode  centerplane.  In  accordance 
with  pinched  beam  diode  theory,  the  excess  negative  space  charge  partially 
suppressed  electron  emission  while  enhancing  ion  emission  to  the  levels,  I  ■ 
0.75  MA  and  I^  ■  0.50  MA.  This  gives  an  ion-to-electron  current  ratio  of  0.67 
compared  to  a  theoretical  prediction  of  0.7  arrived  at  via  Eq.  (4).  The 
overall  agreement  is  excellent.  This  encouraged  one  further  step  of 
investigation.  It  was  found  to  be  true  with  low  inpedance  axial  diodes  that 
the  substitution  of  an  electron-transparent  anode  foil  in  place  of  the  solid 
anode  would  significantly  boost  the  diode's  ion  production  efficiency.  That 
same  substitution  was  tested  here  and,  indeed,  complementary  results  were 
observed.^  As  illustrated  in  Fig.  6,  a  dense  sheath  of  reflexing  electrons 
formed  along  most  of  the  anode  foil  surface.  This  additional  negative  space- 
charge  enhanced  ion  emission  in  the  half-diode  to  0.7  MA  while  further 
reducing  the  electron  current  to  0.6  MA.  The  new  ratio,  I^/Ie,  of  1.17  was 
now,  as  expected,  far  in  excess  of  the  solid  anode  prediction.  As 
anticipated,  the  correlation  of  detailed  diode  physics  between  the  radial  and 
axial  geometries  was  very  strong.  Further  proof  of  this  can  be  found  in  the 
sample  electron  trajectory  plots  of  Fig.  7.  In  order  to  generate  these  plots, 
the  electric  and  magnetic  field  structures  self-consistent  with  the  steady- 
state  particle  flows  of  the  diode  were  "frozen"  over  the  computational  mesh. 
Into  this  field  structure,  electron  streams  were  injected  at  sample  electron 
emission  points.  The  plots  show  the  resultant  stream  trajectories  for  five 
sample  cathode  emission  points.  Their  behavior  compared  closely  to  that  of 
their  counterparts  in  axial  geometry.^  These  simulations  had  therefore  shown 


that  there  were  no  unexpected  differences  between  actual  electron  flow  in  the 
radial  diode  configuration  and  that  postulated  in  the  decoupling  theory.  The 
overall  numerical  results  for  this  radial  diode  are  summarized  in  Table  1. 


Encouraged  by  these  results,  the  next  step  taken  was  the  initiation  of  an 
experimental  test  of  the  decoupling  mechanism  on  the  Gamble  II  pulsed  power 
generator. 


Table  1:  SHEPERD  Diode  Results 


Anode  Type 

Ion  emission 

Ie(MA) 

If  (MA) 

ni 

Z(n)(full  Diode) 

Solid 

None 

0.85 

0 

0 

0.88 

Solid 

Partial 

1.075 

0.075 

0.07 

0.65 

Solid 

Partial 

1.03 

0.14 

0.12 

0.64 

Solid 

Full 

0.75 

0.50 

0.40 

0.60 

Foil 

Full 

0.61 

0.70 

0.53 

0.57 

lO 

The  details  of  the  experimental  setup  may  be  found  elsewhere.  Only  the 

immediately  relevant  points  will  be  related  here.  To  begin  with,  the  nature 

of  the  Gamble  II  machine  argued  against  a  test  of  the  full  radial  diode 

geometry.  To  ensure  symmetric  power  flow  into  both  halves  of  such  a  diode,  a 
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tri-plate  vacuum,  power  transmission  line  should  be  employed.  This  can  best 
be  accomplished  at  Sandia  National  Laboratories  on  their  PBFA-I  device. 

Gamble  II,  on  the  other  hand,  embodies  a  conventional,  coaxial  power  feed  to 
its  diode  load.  Barring  the  use  of  complicated  power  feed  convo lutes,  and 
wishing  to  capitalize  on  the  power  characteristics  of  the  Gamble  II  machine, 
it  was  decided  to  construct  and  test  a  diode  which  physically  represented  only 
half  of  the  full  radial  pinch-reflex  configuration.  The  diode  geometry  was 


spile  at  its  centerplane  in  exact  analogy,  by  coincidence,  to  the  numerical 
simulation  model  but,  of  course,  with  a  perfectly  absorbing  "centerplane", 
"right"  boundary  instead  of  a  perfectly  reflecting  one.  This  experimental 
configuration  is  depicted  in  Fig.  8.  Hie  power  feed  is  from  the  left.  The 
anode  and  cathode  structures  are  clearly  shown  to  approximate  scale.  In  this 
geometry,  the  diode  radius  remained  limited  to  between  6.4  and  7.0  cm  while 
the  A-K  gap  and  the  all-important  half-width,  L,  could  be  varied  between  0  to 
1.0  cm  and  5,  10  and  15  cm  respectively.  Two  sets  of  preliminary  findings 
were  derived.  First,  with  the  A-K  gap,  D,  fixed  at  4.75  mm,  electrode 
structures  of  three  different  lengths,  L  ■  5,  10,  15  cm  were  tested.  As  shown 
in  Fig.  9a,  the  impedance  was  found  to  be  quite  Insensitive  to  the  ion- 
efficiency  -determining  parameter,  L.  On  the  other  hand,  when  a  fixed  L  of  10 
cm  was  tested  with  gap  widths,  D,  of  3.7,  4.8,  and  5.9  mm,  the  impedance  still 
scaled  linearly  with  D  as  expected  from  the  theory  (see  Fig.  9b).  Thus,  for  a 
predetermined,  machine  dependent  impedance,  only  the  radius  and  A-K  gap  of  a 
diode  need  remain  fixed  while  the  experimenter  is  free  to  vary  L  to  achieve  a 
full  range  of  I^/Ig  ratios  as  desired.  For  this  particular  configuration, 
optimization  of  light  ion  production  was  the  primary  consideration.  Thus, 
anode  half-widths  of  up  to  three  times  the  diode  radius  were  of  interest  to 
boost  the  ion  current  fraction  to  about  three  times  that  normally  realizable 
in  an  axial  diode  of  the  same  impedance.  For  other  applications,  the  converse 
problem  of  minimizing  ion  production  is  the  major  objective.  The  alternate 
solution  is  described  in  the  following  section. 

III.  ENHANCED  ELECTRON  CURRENT 

A.  Radial  Geometry 

To  achieve  a  minimization  of  relative  ion  current  in  a  diode  of  a  certain 


impedance,  the  same  radial  diode  concept  as  that  shown  in  Fig.  2  and  discussed 
for  the  case  of  ion  enhancement  in  Section  II  can  be  employed.  First,  the 
diode  current  and,  hence,  its  impedance  is  fixed  to  some  desired  value  by 
adjusting  the  R/D  ratio  found  in  Eq.  (1).  Then  the  ratio  of  1^  to  Ig  is 
minimized  by  shrinking  the  half-width,  L,  using  Eq.  (4)  with  L  substituted  for 
R.  For  a  certain  set  of  experimental  parameters  that  were  of  interest  for  a 
specific  appliction,  the  diode  geometry  depicted  in  Fig.  10  was  chosen  for 
study.  Note  the  similarity  to  the  ion  diode  of  Fig.  3.  Since  a  focused  flow 
of  electrons  is  now  the  prime  product  of  interest,  the  positions  of  the 
cathode  and  anode  have  been  exchanged.  To  attain  a  low  impedance  (i.e.,  a 
high  current  flux),  a  cathode  radius  of  10.0  cm  was  used  in  conjunction  with  a 
relatively  small  A-K  gap  of  0.2  cm.  This  sets  R/D  to  a  value  of  50.  At  the 
same  time,  L  is  kept  small,  0.6  cm,  leaving  an  L/D  of  only  3.  Simulations 
using  the  PREMAS  code  were  again  initiated,  as  for  the  previous  diode,  seeking 
equilibrium  field  structures  and  particle  flow  characteristics.  For  this 
diode,  however,  the  voltage  was  fixed  at  1.0  MV. 

To  check  the  details  of  diode  performance,  two  separate  situations  were 
modeled.  In  the  first,  only  electrons  were  permitted  in  the  system.  In  the 
second,  both  electron  and  ion  emission  were  allowed.  The  numerical 
computational  region  used  for  both  situations  was  very  similar  to  that  used 
for  the  ion  diode  of  Section  II.  This  new  region  is  depicted  in  Fig.  11.  It 
extended  radially  from  the  solid,  cylindrical  anode  surface  at  R  ■  10.0  cm  out 
to  R  ■  10.5  cm,  and  axially  from  the  left,  solid,  planar  anode  surface 
designated  as  Z  ■  0.0  over  to  the  diode  centerplane  at  Z  ■  0.9  cm.  As  shown 
in  the  figure,  it  was  spanned  radially  by  50  cells,  each  0.01  cm  wide,  and 
axially  by  64  cells  which  were  0.0125  cm  wide.  This  gave  a  total  cell  count 
of  a  modest  (64+2)  x  (50+2)*  3432  cells.  Electron  emission  was  permitted  from 


all  three  surfaces  of  the  protruding,  0.2  cm  wide,  0.3  cm  long  cathode  ring  as 
well  as  from  the  solid  cathode  surface  radially  recessed  0.3  cm  from  the  face 
of  the  ring.  In  the  second  simulation,  the  emission  of  protons  was  "turned 
on"  over  the  full  0.9  cm  length  of  the  cylindrical  anode  face. 

In  running  the  "electrons-only”  case  It  was  again  of  interest  to  see  how 
closely  the  current  flow  mimicked  that  normally  observed  in  conventional 
axial  diodes.  At  the  start  of  the  simulation,  the  equipotential  lines  in  the 
source-free  diode  gap  were  plotted  (see  Fig.  12a).  It  is  clear  from  this 
figure  that  2-0  enhancement  of  the  normal  electric  field  at  the  edges  of  the 
cathode  face  is  present  and  that,  therefore,  enhanced  electron  emission  can  be 
expected  there.  After  electron  emission  is  initiated  and  steady  state 
operation  achieved,  the  equipotentials  readjust  themselves  to  accommodate  the 
accumulations  of  negative  space  charge  (see  Fig.  12b).  The  Gaussian  emission 
of  electrons  at  the  cathode  surface  forces  the  normal  electric  field  there  to 
zero,  effectively  "pushing  back"  the  equipotentials  in  that  vicinity.  At  the 
same  time,  the  contours  are  conpressed  near  the  opposing  anode  surfaces. 

To  fully  characterize  the  flow  of  electrons  in  the  A-K  gap,  four  other 
specific  diagnostics  were  output.  First  of  all,  the  electron  current  density 
hitting  the  anode  cylinder  as  a  function  of  Z  is  plotted  in  Fig.  13a.  The 
scatter  of  data  points  shows  a  clear  choppiness  to  the  current  flow  but  also 
indicates  a  sharp  peak  at  about  Z  ■  0.35  cm,  a  rough,  high  plateau  between  Z  » 
0.45  to  Z  -  0.60  cm  and  a  lower  peak  near  the  diode  centerplane  between  Z  * 
0.85  to  Z  -  0.90  cm.  Integration  over  the  half-diode  yields  an  electron 
current  of  1.51  MA  which  is  significantly  but  not  dramatically  higher  than  the 
theoretical  critical  current  of  1.18  MA.  Further  insight  into  the  nature  of 
the  electron  flow  is  given  by  the  sample,  equilibrium  electron  orbits  shown  in 
Fig.  13b.  To  generate  this  plot,  the  steady-state  electric  and  magnetic 


fields  were  "frozen"  onto  Che  computational  mesh.  Then  streams  of  electrons 
were  injected  into  the  system  from  preselected,  sample  emission  points  along 
the  cathode  surface.  The  temporal  spacing  between  electrons  in  each  stream 
was  fixed  thereby  permitting  an  inspection  not  only  of  electron  positions  but 
also  of  relative  velocities.  (Hear  from  the  picture  is  the  effective 
"focusing"  of  the  electron  flow  from  the  outer  cathode  surface  down  to  the 
cylindrical  anode  between  Z  =  0.3  and  Z  *  0.4  cm  accounting  for  the  peak  in 
Jg(Z)  there.  The  flow  from  the  cylindrical  face  of  the  protruding  cathode 
ring  likewise  hits  the  anode  between  Z  ■  0.4  and  Z  -  0.6  cm,  explaining  the 
roughly  constant  current  density  plateau  there.  Finally,  a  strong  confluence 
of  the  sample  electron  streams  takes  place  near  the  cathode  centerplane.  The 
smooth,  uniform  nature  of  the  magnetic  bending  and  focusing  that  the  streams 
experience  appears  qualitatively  obvious  from  Fig.  13b.  A  quantitative 
assessment  of  this  bending  is  given  by  Fig.  13c.  The  tangent  of  the  angle 
toward  the  centerplane  made  by  an  average  electron  momentum  vector  with  the 
surface  normal  vector  was  plotted  for  each  mesh  point  as  a  function  of  the 
axial  dimension,  Z.  The  average  tangent  values  are  seen  to  vary  only  slightly 
around  1.0  over  almost  the  entire  anode  inpact  surface.  This  indicates  that 
the  average  electron  at  almost  every  value  of  Z  hit  at  about  45°.  Plotting  an 
overall  angular  distribution  function  of  all  electrons  hitting  the  anode  over 
a  prolonged  period  of  time  yields  Fig.  13d.  This  shows  that  very  nearly  every 
electron  hit  within  a  relatively  small  "window"  about  45°.  It  may  be  presumed 
that  the  predominant  cause  for  the  near-zero  average  tangent  values  between  Z 
■0.8  and  Z  ■  0.9  cm  was  the  crossover  of  electrons  from  the  mirror  half -diode 
on  the  other  side  of  the  centerplane.  Electrons  from  the  other  side  hitting 
at  -45°  canceled  the  tangent  values  of  those  hitting  at  4-45°. 

At  this  point,  ions  (specifically,  protons)  were  introduced  into  the 


simulation.  This  may  occur  physically  if  a  hydrogen-rich  plasma  is  created 
along  the  anode  surface.  Numerically,  icn  emission  is  achieved  quite  simply 
by  monitoring  the  normal  component  of  the  electric  field  along  a  predetermined 
portion  of  the  anode  surface  and  by  introducing  sufficient  positive  charge 
carriers  along  that  surface  to  zero  that  field  component.  For  the  purposes  of 
this  simulation,  such  ion  emission  was  "turned  on"  from  Z»0.1  cm  all  the  way 
through  the  centerplane  of  the  diode.  Hie  results  were  dramatic  and  closely 
paralleled  what  takes  place  in  a  standard  axial  pinched  beam  diode.  Consider 
the  situation  after  equilibrium  was  reached  for  this  combined  electron-ion 
case. 

First  of  all,  compare  the  new  steady-state  electrostatic  potential 
contour  plot  shown  in  Fig.  14a  with  its  electrons-only  counterpart  in  Fig. 

12b.  Note  that  the  contours  that  had  hugged  the  anode  suface  in  the  earlier 
plot  are  now  "pushed  back"  by  the  emitted  ion  space  charge.  The  electric 
field  is  zeroed  along  all  emission  surfaces  resulting  in  an  intensification  of 
the  radial  field  midway  in  the  radial  A-K  gap.  Furthermore,  the  potentials 
along  the  cathode's  axial  gap  face  have  also  been  swept  further  back 
indicating  an  intensification  of  the  electron  current  emission  there.  In 
fact,  as  can  be  expected  due  to  the  injection  of  massive  amounts  of  positive 
space  charge,  the  overall  diode  electron  current  has  increased  by  over  30%. 

The  net  average  electron  current  reached  an  equilibrium  value  of  4.19  MA  in 
response  to  the  0.31  MA  of  ion  current  (or  2.1  and  0.16  for  the  half-diode). 

An  even  more  dramatic  change  is  observed  in  the  axial  profile  of  electron 
current  density  impacting  the  cylindrical  anode  face.  This  plot  is  shown  in 
Fig.  14b.  It  bears  no  resemblance  to  the  corresponding  plot  for  electrons- 
only  flow  in  Fig.  13a.  However,  the  sharply  peaked  nature  of  this  profile  is 
precisely  what  one  could  expect  from  classic  pinched-beam  diode  operation. 


The  electron  current  density  hitting  the  anode  peaks  at  the  center  due  to  the 
radially  pinched  flow.  The  central  peak  of  over  500  kA/cm^  is  more  than  50 
times  more  intense  than  Je  opposite  the  cathode  tip  (whose  axial  position  is 
indicated  by  two  dashed  lines  at  Z“0.3  and  0.5  cm).  Integrated  over  the 
entire  anode,  the  net  electron  current  is  4.19  MA,  compared  with  3.02  MA  for 
the  electrons-only  case. 

If,  as  in  the  electrons-only  case,  most  of  the  electrons  were  streaming 
from  the  cathode  to  the  anode  in  simple,  straight  arcs,  then  one  might  now 
expect  a  peak  in  the  electron  angular  impact  distribution  at  some  large 
positive  value.  Specifically,  since  most  electrons  will  originate  at  the 
cathode  ring  tip  (z  ~  0.3  cm)  and  will  impact  the  anode  at  about  z  *  0.9  cm, 
the  mean  tangent  could  be  predicted  to  be  about  (0.9  -  0.3)/0.2  or  +3.  In 
fact,  this  is  not  the  case.  The  actual  angular  distribution  shown  in  Figure 
14c  depicts  a  gentle  peak  well-centered  around  tan  9  *  0.0.  In  search  of  a 
reason  for  such  a  distribution,  one  first  may  examine  the  axial  profile  of 
mean  impact  tangents  plotted  in  Figure  14d.  Below  z  -  0.7  cm,  electrons  end 
their  cross-gap  trajectories  at  consistently  sharp  angles.  Near  the 
centerplane,  above  z  »  0.75  cm,  however,  the  mean  tangents  are  clustered 
around  zero  as  they  were  in  Figure  13c  for  the  electrons-only  case.  Given  the 
sharp  current  density  peak  around  the  center,  it  is  only  reasonable  that  the 
overall  mean  tangent  also  be  about  zero.  Then  the  question  becomes,  why  are 
the  mean  tangents  at  the  centerplane  all  about  zero?  The  simple  argument 
regarding  electrons  crossing  over  from  the  "other  side"  of  the  centerplane  and 
cancelling  large  positive  tangents  with  their  negative  tangents  does  not  apply 
here.  If  that  were  the  case  here,  then  the  distribution  may  indeed  have  an 
average  of  zero  but  it  would  be  approximately,  symmetrically  double-peaked, 
perhaps  with  one  peak  around  +3.0  and  another  near  -  3.0.  The  true  answer 


regarding  Che  peak  at  zero  appears,  instead,  Co  be  Che  electron  orbit 
randomization  illustrated  in  Figure  15.  In  this  snapshot  of  steady-state 
trajectories,  the  scrambling  effect  of  E  x  B  drift  orbits  for  this  pinched 
beam  situation  is  manifest.  Note  the  multiple  figure-eight  paths  about  the 
centerplane  which  numerous  electrons  execute  before  ending  their  maneuvers  at 
the  anode  surface.  Such  orbit  behavior  is  also  observed  for  axial  pinched- 
beam  diodes.  This  orbit  randomization  fully  justifies  a  zero  mean  for 
tan  9  .  In  fact,  the  explanation  is  so  strong  that  one  mist  ask  why  the 
distribution  of  Figure  14c  is  not  precisely  centered  at  0.0.  The  reason  for 
that  asymmetry  is  simply  that  Figure  14c  represents  a  time  integration  of 
electrons  impacting  over  just  half  of  the  axial  anode  surface.  If  this 
integration  had  been  extended  over  infinite  time  and  over  the  entire  anode 
surface,  the  angular  distribution  would  be  absolutely  symmetric  around  zero. 

Table  2:  Enhanced  Electron  Radial  Diode  Results 

Species  Present  Ie(MA)  I^MA)  Z(G)(Full  Diode) 

Electrons  Only  1.51  0  0  0.33 

Ion  and  Electrons  2.08  0.16  0.07  0.22 


The  results  of  these  two  simulation  runs  are  summarized  in  Table  2.  Note 
that  the  currents  presented  in  the  table  are  for  the  half-diode.  Multiply  by 
two  to  obtain  the  full  diode  results.  Of  primary  interest  is  the  boosting  of 
diode  electron  current  from  3  MA  (*1.25  Icrit)  to  4.15  MA  (*1.7  Icrit)  through 


the  introduction  of  ion  emission  in  the  radial  A-K  gap.  This  enhancement  of 
electron  current  is  highly  desirable  for  specific  applications  such  as  x-ray 
production  through  bremsstrahlung.  What  is  not  desired  is  a  "loss"  of  too 


much  diode  power  to  unnecessary  ion  production.  Minimization  of  ion  current 
motivated  this  particular  choice  of  diode  geometry.  The  ratio,  R/D,  was  set 
at  10.0/0.2  ■  50  to  boost  Icr^t  (and  thereby,  Ig)  via  Eq.  (3).  At  the  same 
time,  the  diode  gap  aspect  ratio,  L/D,  was  set  to  0.6/0. 2  3  3  so  that  Eq.  (4) 
(with  L/D  substituted  for  R/D)  yields  a  species  current  ratio  of 
If/ I  2  0.07.  The  simulation  found  a  ratio  of  0.078  which  amounts  to 
excellent  agreement  with  theory.  Thus,  the  "cost"  of  boosting  Ie  from  1.25  to 
1.7  times  I  may  be  considered  “cheap"  in  terms  of  wasted  ion  current. 

There  is  a  problem,  however,  with  this  configuration.  The  radial 

geometry  results  in  a  radial  flow  of  electron  current  with  a  resultant  radial 

flux  of  the  bremsstrahlung  x-rays.  Thus  the  object  to  be  exposed  to  the  x- 

rays  (say,  a  patient  for  medical  applications)  must  be  inserted  into  the 

hollow  anode  volume.  This  places  an  obvious  upper  bound  on  the  physical  size 

of  the  object  to  be  irradiated.  In  addition,  there  are  high  voltage  hazard 

considerations  to  be  dealt  with  when  placing  anything  into  the  diode  region. 

These  difficulties  could  be  alleviated  if  this  high  fluence  diode  could  be 

configured  in  an  axial  rather  than  radial  geometry.  Such  a  change  wou'.i  allow 

a  predominant  x-ray  flux  away  rrom  the  diode  structure  which  could  be  focused 

or  defocused  to  irradiate  a  small  or  an  extended  target.  Clearly,  a 

conventional  axial  diode  design  will  not  achieve  this  objective.  The 

solution,  a  magnetically  insulated  splitter  (or  MIS)  diode,  was  suggested  by 
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researchers  at  Maxwell  Laboratories.  The  essence  of  this  device  relied  on  a 
rerouting  of  the  diode  return  current  through  the  inside  of  the  hollow 
cathode,  thus  altering  the  strength  of  the  diode  B  felt  by  the  electrons  as 
they  transit  the  A-K  gap.  The  MIS  concept  and  a  simulation  of  the  basic 
device  are  presented  in  the  following  section. 


B.  Axial  Geomet 


The  geometrical  translation  of  the  enhanced  electron  current  radial  diode 
discussed  in  the  previous  section  to  a  practical  axial  diode  conf iguaration 
was  straightforward.  The  key  element  which  spelled  success  for  the  radial 
diode  was  the  physical  positioning  of  the  imaginary  B  ■  0  surface  within  the 
A-K  gap  in  such  a  way  that  it  approximately  separated  one  half  of  the  net 
diode  electron  current  from  the  other  half.  Since  the  B  *  0  line  traced  on 
the  anode  surface  represents  the  focus  of  the  pinched  electron  flow,  its 
position  determines  the  approximate  electron  "dwell  time"  in  the  A-K  gap  and 
thus  impacts  on  the  ion-to-electron  current  ratio  in  accordance  with  Eq. 

(4).  To  achieve  this  effect  in  axial  geometry,  the  radial  geometry  is 
schematically  "rotated".  The  radial  axis  is  thus  transposed  with  the  axial 
axis.  The  physical  equivalence  of  the  two  orientations  is  achieved  by  routing 
half  of  the  net  diode  current  through  a  conductor  located  at  a  radius  inside 
the  inner  radius  of  the  hollow  cathode.  This  is  accomplished  in  the  MIS 
concept  device  previously  mentioned.  In  this  new  geometry,  R/D  is  once  again 
replaced  by  L/D  in  the  determination  of  I^/Ie.  However,  L  is  now  defined  as 
simply  the  halfwidth  of  the  hollow  cathode  cylinder  (i.e.,  one  half  of  the 
distance  between  the  cathode's  inner  and  outer  radius).  Once  again,  a  cathode 
with  small  L  can  be  used  to  restrict  I^/Ie  while  simultaneously  keeping  R 
large  with  respect  to  D  in  order  to  preserve  a  high  total  electron  current. 

To  test  this  concept  for  axial  diode  electron  enhancement,  the  geometry 
depicted  in  Fig.  16  was  modeled.  This  cathode  was  taken  as  a  solid  ring  with 
a  mean  radius  of  13.0  cm  and  with  axially  protruding  inner  and  outer  radial 
edges.  These  edges  protrude  about  one-eight  inch  (~  0.34  cm)  out  from  the 
solid  ring.  This  ring  cathode  is  almost  completely  surrounded  by  a  hollow 
toroidal  anode  shell.  This  shell  is  radially  centered  on  the  ring  so  that  its 


mean  radius  is  also  13.0  cm.  A  radial  A-K  gap  of  3.0  cm  is  maintained  on 
either  side  of  the  ring.  Since  the  ring  has  a  radial  thickness  of  1.6  cm, 
this  implies  a  radial  anode  thickness  of  7.6  cm.  Finally,  the  effective  A-K 
gap  (after  plasma  gap  closure)  is  set  at  0.3  cm.  In  order  to  mimick  an 
experimentally  interesting  device,  the  constant  diode  voltage  is  set  to  1.5 
MV.  In  constructing  the  specific  numerical  model,  only  the  volume  spanning 
the  full  radial  gaps  and  frontal,  axial  gap  surrounding  the  two  cathode  tips 
is  of  interest.  Thus,  the  rear  planar  boundary  of  the  computational  region  is 
aligned  with  the  recessed  back  surface  in  the  interior  of  the  cathode  ring. 

The  opposite  axial  planar  boundary  as  well  as  the  lower  and  upper  radial 
boundaries  are  naturally  made  to  coincide  with  the  inner  surfaces  of  the 

I 

surrounding  anode  shell.  The  R.-Z  cross  section  of  the  computational  region 
thus  measures  0.64  by  7.59  cm. 

The  resulting  simulation  configuration  is  shown  in  Fig.  17.  Numerical 

1 

cell  spacings  of  AR  x  AZ  *  0.03  x  0.02  cm  are  chosen  to  fill  the  modeling 
region.  This  results  in  32  equispaced  cells  axially  and  253  radially. 

Allowing  a  monolayer  of  "guard  cells”  surrounding  the  entire  perimeter  of  the 
region  gives  a  total  cell  count  of  34x255*8670.  The  protruding  cathode  tips 
are  each  represented  by  conducting  "surfaces"  seven  cells  wide  radially  and  17 
cells  axially.  They  are  maintained  at  the  cathode  potential  of  zero  as  is 
that  segment  of  the  left  boundary  which  is  sandwiched  between  them.  The 
remaining  two  segments  of  the  left  boundary  separating  the  anode  and  cathode 
have  fixed  logarithmic  gradings  of  potential  between  zero  and  4-1.5  W.  When 
the  simulation  is  underway,  electron  emission  is  permitted  only  from  the  two 
seven-cell-wide  faces  of  the  protruding  cathode  tips  and  from  the  inner 
surface  of  the  cathode  recess.  Two  separate  electrons -only  steady  states  are 
simulated.  In  the  first,  no  return  current  is  allowed  to  flow  Inside  the 
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cathode  radius.  In  the  second,  exactly  half  of  the  diode  current  is  directed 
along  the  lower  axial  boundary  which  represents  the  inner  conducting  surface 
of  the  anode  shell.  In  the  final  run,  the  50-50  division  of  the  diode  return 
current  is  maintained  while  proton  emission  is  "turned  on"  along  the  circular 
strip  of  anode  surface  opposite  the  ring  cathode  face. 

The  simulation  region  is  initialized  at  t=Q  with  no  charges  present  and 

only  a  fraction  of  the  peak  diode  voltage  imposed.  The  electrostatic 

equipotential  contours  which  result  in  this  configuration  are  shown  at  10% 

intervals  in  Fig.  18.  Note  the  exceptionally  tight  compression  of  the 

contours  in  the  vicinity  of  the  two  tips  with  the  resultant  intensification  of 

normal  electric  field  values  there.  This  is  a  manifestation  of  the  two- 

dimensional  electric  field  enhancement  typical  in  such  hollow  cathode 

devices.  As  expected,  as  soon  as  electron  emission  is  "turned  on"  intense 

streams  of  electron  current  are  ejected  from  the  two  protruding  tips.  At  the 

same  time,  the  low  electric  fields  in  the  interior  region  between  the  two  tips 

cause  more  than  an  order-of -magnitude  lower  emission  from  the  cathode  ring's 

recessed  back  face.  This  relatively  slow  moving  mass  of  electron  space  charge 

in  this  interior  region  sweeps  back  the  electrostatic  potential  contours  there 

resulting  in  a  further  decrease  of  the  already  low  electric  field  and  the 

consequent  reduction  of  electron  emission  there  to  a  mere  trickle.  This 

field-modification  process  is  clearly  visable  in  the  equipotential  contour 

plot  shown  in  Fig.  19a  as  compared  with  that  of  Fig.  18.  In  this  new  figure, 

steady-state  electron  flow  is  fully  established.  One  thousand  timesteps  of 
-12 

0.5x10  second  each  have  passed  and  this  translates  to  over  five  typical 
electron  A-K  gap  crossing  times.  The  regions  of  high  concentrations  of 
electron  space  charge  manifest  themselves  not  only  in  the  just -mentioned 
sweep-out  of  contours  from  the  cathode  ring  interior  but  also  in  the  double- 
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horned  distortion  of  all  of  the  contours  where  the  two  principle  electron 
current  streams  are  localized.  Keep  in  mind  that  all  of  the  diode  return 
current  flows  along  the  outside  (upper)  radial  boundary  in  this  run. 

Therfore,  except  for  effects  associated  with  the  protrusion  of  the  anode 
surface  inside  the  inner  radius  of  the  cathode  ring,  the  behavior  of  this 
device  should  closely  resemble  that  of  a  standard,  axial,  hollow  cathode 
diode.  It  does. 

Just  as  in  the  previous  set  of  simulations  for  the  radial  enhanced  I 

diode,  so  too  here  detailed  diagnostics  were  performed  at  the  anode  plane  in 

order  to  obtain  a  comprehensive  characterization  of  the  electron  beam.  The 

first  diagnostic  of  Fig.  19b  shows  the  radial  profile  of  the  electron  current 

density  striking  the  anode.  As  a  convenient  reference,  the  radial  positions 

of  the  two  protruding  cathode  tips  on  the  opposite  side  of  the  diode  are 

indicated  on  the  plot  with  dotted  lines.  The  lines  clarify  the  degree  of 

radial  deflection  the  beams  have  undergone  in  their  traverse  of  the  axial  a-K 

gap.  In  a  standard  hollow  cathode  diode  with  suppressed  io'*  flow,  would 

expect  the  hollow  electron  beam  to  have  been  pinched  radially  inward  by  a 

distance  about  equal  to  the  A-K  gap.  This  is  consistent  with  the  critical 

current  model  explained  in  Section  I.  A  typical  electron  may  be  pictured  to 

have  executed  a  90°  circular  arc  in  the  gap  transit.  Such  deflections  appear 

in  Fig.  19b.  Both  beam  components  are  pinched  about  0.3  cm  inward  toward 

their  common  centerline  after  their  0.3  cm  axial  traverse.  Each  peak  emitted 

by  its  respective  cathode  tip  is  sharply  pronounced  with  the  inner  beam 

carrying  more  current  than  the  outer  beam.  The  outer  cathode  tip  emits  less 

electron  current  since  it  is  immersed  in  the  Ba  field  generated  by  the  inner 

cathode  tip  and  beam.  This  strong  B  acts  to  partially  insulate  emission  on 

0 

9  9 

the  outer  tip.  As  an  estimate  for  the  expected  net  current,  one  may  use  Eq. 
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(3)  with  Sg  ~  1,  Ye  *  3.935,  R  ■  13.0  cm,  and  D  =  0.3  cm.  This  yields  lcr^t  * 
1.45  MA.  If  the  two  costreaming  hollow  electron  beams  were  indeed 
independent,  -.his  would  predict  a  net  diode  current  of  about  2  Icr^t  *  2.9 
MA.  In  fact,  as  just  seen,  the  outer  beam  is  magnetically  coupled  to  the 
inner  one  and  experiences  a  current  reduction  because  of  it.  This  helps  to 
explain  the  net  observed  diode  electron  current  of  2.16  MA  for  this  case  in 
steady  state. 

The  other  major  beam  diagnostics  are  concerned  with  the  local  and  global 
angular  distribution  of  the  electrons  incident  on  the  anode  face.  First,  a 
radial  scan  of  impact  angles  as  a  function  of  the  anode  radius,  r,  yields  the 
plot  of  Fig.  19c  At  each  0.03  cm  wide  data  cell,  the  average  is  taken  of  the 
tangent  of  their  impact  angle,  0.  As  expected  from  the  radial  pinching,  all 
of  the  angles  are  positive.  Furthermore,  if  all  of  these  data  points  were 
equally  weighted,  one  would  calculate  an  overall  mean  inpact  tangent  of 
slightly  over  2.0.  However,  all  points  are  not  of  equal  weight  since  the 
current  density  profile  is  far  from  flat.  There  are  essentially  two  well- 
defined  hollow  beams.  From  this  plot  one  may  assume  that  the  inner  beam  is 
represented  by  the  cluster  of  data  points  around  tan  0  s  1.0  below  r  *  12.3 
cm.  Thus,  typical  electrons  from  that  beam  are  striking  the  anode  at  about 
45°.  Similarly,  the  trace  of  points  above  tan  0  a  2.0  and  r  *  13.0  cm  are  due 
to  the  electrons  of  the  outer  beam.  They  experience  stronger  magnetic  bending 
because  they  feel  the  Bg  generated  by  the  inner  beam,  thus  they  strike  the 
anode  at  shallower  angles  (larger  0's  with  respect  to  the  surface  normal). 
Weighting  these  mean  tangents  by  their  respective  current  densities  from  Fig. 
19b  yields  the  global  angular  distribution  plot  of  Fig.  19d.  It  is  the 
behavior  of  the  more  intense  inner  beam  that  dominates  the  overall  physics. 

The  distribution  peaks  slightly  below  tan  0  ■  1.0  similarly  to  that  observed 


for  the  corresponding  radial  diode  case. 

At  this  point  in  the  simulation  (t  =  1000)  exactly  half  of  the  net  diode 
current  at  any  given  time  is  rerouted  along  the  inner  (r  *  9.21  cm)  surface  of 
the  anode  "can".  The  effect  on  the  dual-beam  electron  flow  is  exactly  as 
expected.  By  t  =*  2000,  a  near  equilibrium  for  the  diode  with  the  50-50  return 
current  split  had  been  established.  The  new  configuration  for  the 
electrostatic  equipotential  contour  surfaces  is  shown  in  Fig.  20a.  Even  from 
this  plot  alone  the  new  relative  balance  between  the  inner  and  outer  beams  is 
apparent.  The  former  pinching  toward  the  common  centerline  at  r  *  0  has  given 
way  to  complementary  pinching  toward  the  imaginary  centercylinder  surface 
midway  between  them  at  r  *  13.0  cm.  This  phenomenon  is  confirmed  by  the 
electron  current  density  profile  shown  in  Fig.  20b.  Not  a  single  data  point 
falls  outside  the  inner  and  outer  cathode  ring  radii.  The  two  peak  densities 
from  the  two  respective  cathode  tips  are  still  distinct  and  still  unequal 
although  not  nearly  to  as  great  a  degree  as  in  the  previous  case.  The  lack  of 
perfect  equality  is  traced  simply  to  the  inequality  of  the  tip  radii.  The 
inner  tip  is  at  r  *  12.2-12.4  cm  while  the  outer  one  is  at  r  *  13.6-13.8  cm. 
The  currents,  though  magnetically  decoupled,  cannot  be  equal. 

An  even  more  dramatic  and  pleasing  manifestation  of  this  new  centered 
pinching  mechanism  between  the  cathode  tips  appears  in  Fig.  20c  of  the  radial 
profile  of  mean  impact  tangents  along  the  anode  face.  The  profile  is  nearly 
perfectly  balanced  with  opposing  impacts  on  opposite  sides  of  r  *  13.0  cm. 
Electrons  from  the  inner  tip  impact  at  negative  angles  (tan  9  ~  -1)  while 
those  from  the  outer  tip  strike  at  positive  angles  (tan  9  ~  +1).  In  deriving 
an  overall  distribution  from  this  radial  angular  profile  once  again  the  radial 
density  profile  oust  provide  the  relative  weighting.  Synthesizing  Figs.  20b 
and  20c  would  lead  one  to  expect  a  global  distribution  function  with  two 
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distinct  peaks;  one  at  tan  9  =  -1  and  one  near  tan  9  *  +1.  That  is  indeed 
what  results  in  Fig.  20d.  In  this  case  too,  the  symmetry  about 
the  tan  9  ■  0.0  is  far  from  perfect.  Again,  this  asymmetry  may  be  attributed 
to  the  unequal  tip  radii.  The  inner  tip  normally  emits  more  current  than  the 
outer  when  all  of  the  return  current  flows  along  the  outer  anode  radius. 
Therefore,  routing  half  of  the  total  diode  current  along  the  inner  radius  will 
not  totally  cancel  the  inner  electron  beam  current.  Some  residual  B.  will 
still  remain  to  insulate  some  of  the  electron  emission  there.  Although  the 
decoupling  of  the  two  electron  flows  was  not  complete,  it  was  still  effective 
enough  to  boost  the  net  diode  I&  to  3.4  MA.  This  compares  well  with  the 
2  lcrit  value  of  2.9  MA  calculated  earlier.  The  remaining  effect  to  test  is 
the  influence  of  ion  flow  on  this  configuration.  Of  particular  importance  is 
the  value  of  I^/Ig  which  will  emerge. 

This  final  simulation  of  the  series  was  accomplished  by  allowing  the 
emission  of  protons  along  the  planar  anode  surface  in  a  ring  of  inner  and 
outer  radius  exactly  equal  to  that  of  the  cathode  ring.  Thus,  ion  emission  is 
"turned  on"  between  r  ■  12.2  and  r  -  13.8  cm.  Eighteen  hundred  additional 
timesteps  are  then  allowed  to  pass  before  this  new  steady  state  is  assessed. 
The  assessment  begins  with  an  examination  of  the  altered  equipotential 
contours  in  Fig.  21a.  The  two  streams  still  seem  distinct  but  some  turbulence 
is  apparent.  The  nature  of  this  stream  mixing  is  clarified  by  the  electron 
current  density  profile  of  Fig.  21b.  Only  one  distinct  density  peak  stands 
out  opposite  the  lower  cathode  tip.  Under  the  influence  of  the  neutralizing 
ion  space  charge,  the  dual  electron  beams  have  pinched  together  as  expected. 
Once  again,  since  the  two  beam  currents  are  not  equal,  the  50%  centerline 
return  current  cannot  center  the  pinch  between  the  two  cathode  tips.  The 
electron  current  for  the  new  steady-state  is  about  3.87  MA,  a  significant 
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increase  over  the  electrons-only  case.  The  net  ion  current  measures  about 
0.19  MA  which  translates  to  a  diode  ion  production  efficiency  of  0.047.  By 
comparison,  taking  Eq.  (4)  and  replacing  R/D  with  L/D  *  0. 8/0.3  =  2.67  yields 

I./I  s  0.075  and  a  predicted  ion  efficiency  of  n.  *  (I. /I  )/(l+I. /I  )  *  0.07. 

1  ®  1  16  16 

The  agreement  is  not  bad  considering  the  crudeness  of  the  approximations  made. 

The  next  question  concerns  the  nature  of  the  new  impact  angle 
distribution  for  the  electrons  striking  the  anode  surface.  Part  of  the  answer 
lies  in  Fig.  21c.  The  plot  of  mean  impact  tangents  as  a  function  of  radius 
shows  that  the  electrons  impacting  near  the  current  density  peak  at  f2.21  cm 
have  an  average  impact  tangent  of  zero.  This  is  precisely  analogous  to  the 
previous  cases  studied.  The  presence  of  ion  space  charge  transforms  the 
simple  electron  arc  trajectories  of  the  ion-free  equilibrium  to 
complicated  £  x  5  drift  orbits  which  randomize  their  final  impact  angles.  The 
scattered  data  points  above  r  -  12.4  cm  with  mean  tangents  of  1.0  and  higher 
represent  only  sparse  current  densities  (i.e.,  relatively  few 
macroelectrons) .  This  is  borne  out  again  by  the  overall  angular  impact 
distribution  plot  shown  in  Fig.  21d.  This  net  distribution  of  all  electrons 
in  the  equilibruim  configuration  is  strongly  peaked  about  tan  8  *  0.0.  Ion- 
induced  pinched  beam  electron  flow  has  been  achieved  although  the  pinching  is 
not  centered  between  the  two  cathode  rings.  Once  again,  the  previously 
described  imbalance  between  the  net  emission  of  the  two  separate  rings  causes 
a  net  displacement  of  the  pinch  radially  inward.  This  shifting  of  the  pinch 
indicates  a  mechanism  for  adjusting  the  mean  electron  impact  radius  not  only 
by  altering  the  percentage  of  total  diode  return  current  routed  along  the 
centerline  but  also  by  controlling  the  radial  extent  of  ion  emission  over  the 
anode  face.  Allowing  ion  emission  down  to  smaller  radii  -  say  10  cm  -  would 
quite  probably  draw  the  pinched  cylindrical  beam's  current  density  peak  even 


further  inward.  This  would  be  an  interesting  topic  for  further  study  but  is 


beyond  the  scope  of  this  report.  The  results  from  these  MIS  diode  simulations 
are  summarized  in  Table  3. 

Table  3:  MIS  Diode  Results 

Species  Present  MIS  Current  Ig  1^  Z(D) 


Electrons  Only 

NO 

2.16 

0 

0 

0.69 

Electrons  Only 

50% 

3.40 

0 

0 

0.44 

Ions  and  Electrons 

50% 

3.87 

0.19 

0.047 

0.37 

IV.  GENERAL  CONCLUSIONS 

*  Standard  diode  theory  links  both  impedance  and  ion  efficiency  to  the 
ratio  of  cathode  radius  to  A-K  gap,  R/D.  This  relationship  between  the  two 
parameters  severely  limited  researchers  in  their  choice  of  operating 
regimes.  If  one  desired  a  high  ion  production  diode,  he  was  restricted  to  one 
ohm  or  sub-ohm  levels  of  operation.  Conversely,  a  project  calling  for  a 
minimization  of  ion  current  in  relationship  to  the  electron  current  dictated  a 
boosting  of  diode  impedance  to  a  multi-ohm  level.  These  restrictions  no 
longer  apply.  Table  4  presents  a  synopsis  of  the  major  results  of  all  three 
sets  of  simulations  related  in  this  report.  The  ion-to-electron  current 
ratios  observed  in  the  runs  are  clearly  tracking  the  cathode  half -width  to  A-K 
gap  ratio,  L/D,  and  not  R/D.  (Note  that  the  figures  presented  in  the  table 
are  all  for  half -diodes  of  each  configuration.) 


Table  4:  Summary  of  Results 


Diode  Type 

Enhanced  Ions 

Enhanced 

Electrons 

Geometry 

Radial 

Radial 

Axial 

Voltage  (MV) 

1.5 

1.0 

1.5 

R,  D,  &  L  (cm) 

5.0,  0.2,  5.0 

10.0,  0.2,  0.6 

13.0,  0.3,  0.8 

R/D  and  L/D 

25  and  25 

50  and  3 

43.3  and  2.67 

I^/Ig  from  R/D 

0.707 

1.15 

1.22 

Ii/Ie  from  L/D 

0.707 

0.069 

0.07  5 

Simulation  I^/Ig 

0.67 

0.074 

0.049 

Theor.  Current  (MA) 

1.46 

1.34 

1.59 

Sim.  Current  (MA) 

1.25 

2.25 

2.04 

Theor.  Impedance  (0) 

1.027 

0.746 

0.943 

Sim.  Impedance  (0) 

1.20 

0.44 

0.735 

It  is  also  clear  from  the  table  that  the  enhanced  ion  production  case  is 
not  a  rigorous  test  of  this  technique  since  R  and  L  are  identical.  This 
suggests  an  interesting  follow-on  effort.  Prior  to  the  research  on  the 
SHEPERD  diode,  significant  time  and  energy  had  been  invested  by  the  NRL 
Intense  Beam  group  in  a  DNA-sponsored  project  to  design  an  efficient  ion  diode 
for  the  AURORA  pulsed  power  generator  at  the  Harry  Diamond  Laboratory ^ 
Working  at  voltages  around  3-5  W  and  impedances  of  up  to  40  ohms,  ion 
efficiencies  of  at  best  5%  to  20%  were  realized.  A  radial  diode  configuration 
embodying  a  large  L/D  but  small  R/D  would  hold  great  promise  for  achieving 
50%+  ion  efficiencies  at  20  ohm  levels.  Unfortunately,  such  a  device  would  be 
strongly  elongated  in  the  axial  dimension,  causing  significant  problems  in  the 


design  of  useful  “targets."  Still,  an  intense  flux  of  relatively  high  voltage 
ions  may  be  of  interest  for  specific  applications.  In  addition,  although  an 
axial  geometry  counterpart  is  impossible  (since  L  would  be  restricted  to  be 
less  than  R)  some  intermediate  geometry  such  as  conical  or  hemispherical  may 
be  workable. 
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